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Stability of black holes in Einstein-charged scalar field theory in a cavity 
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Can a black hole that suffers a superradiant instability evolve towards a ‘hairy’ configuration 
which is stable? We address this question in the context of Einstein-charged scalar field theory. 
First, we describe a family of static black hole solutions which possess charged scalar-field hair 
confined within a mirror-like boundary. Next, we derive a set of equations which govern the linear, 
spherically symmetric perturbations of these hairy solutions. We present numerical evidence which 
suggests that, unlike the vacuum solutions, the (single-node) hairy solutions are stable under linear 
perturbations. Thus, it is plausible that stable hairy black holes represent the end-point of the 
superradiant instability of electrically-charged Reissner-Nordstrom black holes in a cavity; we outline 
ways to explore this hypothesis. 

PACS numbers: 04.40.Nr 04.40.Bw 


I. INTRODUCTION 

The belief that a ‘typical’ galaxy hosts a supermassive 
black hole, of mass M ~ is supported by 

dynamical evidence from nearby galaxies and extrapola¬ 
tion of the black hole mass-velocity dispersion relation 
[1]. Recent surveys suggest that a supermassive black 
hole may possess significant angular momentum [2, 3]. 
Black hole spin is conjectured to power relativistic jets 
in quasars through the Blandford-Znajek process [4], fea¬ 
turing an accretion disk and a force-free magnetosphere 

The Blandford-Znajek process is just one example of 
a Penrose process [6] , in which a black hole may liberate 
energy and angular momentum (and/or charge) whilst 
still increasing its horizon area. Penrose processes are 
consistent with - indeed, encouraged by - the second law 
of black hole mechanics [7] and thus, it would appear, 
the second law of thermodynamics [8, 9]. One intriguing 
example of a Penrose process is superradiance [10], in 
which a low-frequency electromagnetic or gravitational 
wave packet is amplified by a black hole (see [11] for a 
recent review). In the ‘black hole bomb’ scenario [12], 
an exponentially-growing instability is stimulated by re¬ 
flecting a superradiant field back onto the black hole. 
In scenarios with light bosonic fields (e.g. axions [13] 
or massive photons [14]), the instability may develop in 
the gravitationally-bound modes of the field [15-25], and 
thus arise spontaneously [26]. Variations on this theme 
involving only electromagnetic fields and accretion disks 
have also been discussed [27]. 

In this article we address a key question: can a super¬ 
radiant instability, pursued into the non-linear regime, 
lead to a new stable hairy black hole configuration? 

Superradiant instabilities appear to pose a challenge 
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to the ‘no-hair’ (Israel-Carter) conjecture [28-30], which 
asserts that a perturbed black hole should settle back 
into a stationary state, changing only a small number of 
parameters (mass M, angular momentum J, and charge 
Q). The conjecture has been codified in a number of the¬ 
orems in asymptotically-flat scenarios involving massless 
scalar, electromagnetic and gravitational fields under cer¬ 
tain minimal assumptions [31-35]; there also exist stabil¬ 
ity theorems for Kerr spacetime [36-38]. Nevertheless, 
as was recently shown in Refs. [39, 40], there exists an 
asymptotically-flat family of Kerr black holes endowed 
with (complex, massive) scalar-field ‘hair’, which reduce 
to well-known ‘solitonic’ boson star solutions in a well- 
defined limit. Crucially, these new solutions lie beyond 
the scope of the no-hair theorems, as the scalar field is 
only helically-symmetric, rather than being stationary 
and axisymmetric. More precisely, the scalar field is only 
invariant under a single Killing field, which is tangent 
to the null generators of the horizon. Novel solutions 
with a single Killing field were described in Ref. [41]. 
For a succinct summary of the relationship between vari¬ 
ous asymptotically-flat scalar-hairy solutions, the no-hair 
theorems, and the violated assumptions, see Table I in 
Ref. [42]. 

In the small-amplitude regime, superradiance arises for 
charged scalar perturbations of the Kerr-Newman space- 
time which have frequency a satisfying a{a — CTc) < 0- 
The critical frequency CTc is given by 

ac = mUln + q4>H, (1-1) 

with UIh = J/2M'^r+ and (fn = Qlf+ the angular fre¬ 
quency and electric potential of the black hole horizon at 
r = r+ = M + — {J/MY — Q^, and m and q the 

azimuthal mode number and charge of the scalar field, 
respectively. At the critical frequency a = ac, linear per¬ 
turbations are stationary: they do not decay or grow (see 
e.g. Ref. [43]). In the limit of small field amplitude, the 
Kerr-scalar solutions in [39, 40] reduce to a Kerr black 
hole {Q = 0) with a co-rotating dipolar (£ = m = 1, 
where I is the total angular momentum mode number) 
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perturbation in the massive scalar field at the critical 
superradiant frequency, a = a^- 

Analysing the stability of the non-linear Kerr-scalar 
solutions (with J 7 ^ 0 , Q = 0 ) is challenging, principally 
because such solutions are only helically-symmetric (as 
well as being numerically-determined, i.e. not known in 
closed form). Here, we consider a simpler spherically- 
symmetric model, with J = 0, Q 7 ^ 0, of scalar electro¬ 
dynamics coupled to gravity [44]. In this scenario, super¬ 
radiance is associated with charge, rather than angular 
momentum. It was shown by Bekenstein (see Ref. [31], 
Sec. IV) that asymptotically-flat finite-energy configura¬ 
tions with charged scalar-field hair are forbidden. In¬ 
stead, we consider an analogue of the black hole bomb 
scenario of Press and Teukolsky [12], in which the black 
hole is enclosed by a reflecting mirror. 

It was shown in Refs. [45-47] that, in the linear (small- 
amplitude) regime, a charged scalar field with a mirror 
on a Reissner-Nordstrom black hole background (J = 0, 
Q 7 ^ 0 ) suffers exponential growth due to superradi¬ 
ance, provided the mirror is sufficiently far from the hori¬ 
zon. Here we consider the progression into the non-linear 
regime. We present charged-scalar black holes which are 
plausible endpoints for the above charge-superradiant in¬ 
stability, and examine their stability under perturbation. 

The outline of this paper is as follows. In Sec. H we de¬ 
scribe our Einstein-charged scalar field model and briefly 
review the instability of Reissner-Nordstrom black holes 
under spherically symmetric charged scalar field pertur¬ 
bations [45-47]. We also present static, spherically sym¬ 
metric black hole solutions with nontrivial charged scalar 
field hair. The charged scalar field has zeros outside the 
event horizon; the reflecting mirror can be situated at 
any one of these zeros. To see if these hairy black holes 
are plausible endpoints of the charge superradiant insta¬ 
bility, in Sec. HI we investigate their stability under lin¬ 
ear, spherically symmetric perturbations. If the mirror 
is located at the first zero of the charged scalar field, we 
present numerical evidence that at least some of the hairy 
black holes are stable. Our conclusions are presented in 
Sec. IV. 


II. BLACK HOLE SOLUTIONS WITH HAIR 


scalar field $. Tensor indices are lowered and raised with 
the metric gab and its inverse g°‘^, and g denotes the met¬ 
ric determinant. Round and square brackets on indices 
denote symmetrized and anti-symmetrized combinations, 
-^(ab) — 2 ^ba^ and X^ab] — 2 ^ba')- 

By varying (2.1), three equation of motions are ob¬ 
tained 


Gab = SirGTab, 

( 2 . 2 a) 

VaF'^^ = J^ 

( 2 . 2 b) 

DaD^^ = 0 , 

( 2 . 2 c) 


alongside the usual Bianchi identities for the Faraday and 
Riemann tensors, X[aFbc] = 0 = '^[aR^^de]- stress- 
energy tensor is given by Tab = T^^^ + where 



rp^ 

^ ah 


= FacF^ ° — -gabFcdF 

= Dl<^*Db)^-^gab 




and the field current is given by 


(2.3a) 

(2.3b) 


^ - $(£)“$)*]. (2.4) 


The current and stress energy are covariantly conserved, 

VaT“ = 0 = 

The charged scalar field $ and vector potential Aa are 
dehned up to the usual gauge freedom, as Fab and Da^ 
are invariant under the mapping 

Aa ^ Aa+q~\,a, (2.5) 


where y is any scalar field. We will make use of this 
freedom both when we consider static equilibrium solu¬ 
tions and time-dependent, spherically symmetric pertur¬ 
bations. 


B. Linear perturbations in electrovacuum 

In the case 4> = 0, the spherically-symmetric solution 
of the field equations (2.2) is the Reissner-Nordstrom 
black hole spacetime, 


A. The model 


ds^ = gabdx°‘dx*’ = —fuNdt^ + fjipfdr^ + ( 2 . 6 ) 


We consider a fully coupled system consisting of grav¬ 
ity, an electromagnetic field and a massless charged scalar 
field. The action is given by: 


5 = 



R 

IQttG 


\FabF'^'> 


d^x, 

( 2 . 1 ) 


where (in appropriate units). 


Jrn = 1 - 


2M 


(r — r+)(r — r_) 

j,2 


and the element of solid angle is 

(Fd? = dO"^ + sin^ 9 dtp^. 


(2.7) 


( 2 . 8 ) 


where Fab = VoAf, — Vf,Aa is the Faraday tensor and 
Da = ^a—iqAa, with Vq the covariant derivative, Aa the 
electromagnetic vector potential and q the charge of the 


The quantities 


r±=M± - g 2 


(2.9) 
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FIG. 1. The real (left) and imaginary (right) part of the frequency ctq of the fundamental massless scalar perturbation of 
Reissner-Nordstrom spacetime with a reflecting mirror, plotted against the mirror radius with (top row) fixed black hole 
charge Q = 0.9 and (bottom row) fixed scalar field charge q = 0.5 (we use units in which the black hole mass M = 1). The 
right-hand plots show the transition points where the imaginary part of ao changes sign from negative (decaying) to positive 
(growing). 


are, respectively, the radii of the outer (event) and inner 
(Cauchy) horizons. 

One may introduce a small-amplitude scalar field 
and neglect the back-reaction on the electromag¬ 
netic and gravitational fields (as Tab and J“ are 
quadratic in the scalar field amplitude). Let us con¬ 
sider a monochromatic, spherically-symmetric perturba¬ 
tion with frequency a 


$ = 


^ 1 


which, via Eq. (2.2c), satisfies 


(P(j) 

dri 



JrN df RN 
r dr 


</. = 0 , 


( 2 . 10 ) 


( 2 . 11 ) 


where the tortoise coordinate r* is defined by dr^/dr = 
fuN- The scalar field perturbation should be ingoing 
at the horizon, that is, regular in a (future) horizon- 
penetrating coordinate system, which implies that 


(j) e as r* —>• —oo. (2.12) 


Imposing a ‘mirror’ boundary condition ipir^) = 0 at 
r = Tm leads to a discrete spectrum of states with, in 


general, complex frequencies (T„. A positive (negative) 
imaginary component of frequency corresponds to expo¬ 
nential growth (decay). The states are labelled with n, 
the number of nodes they possess in the region r < 

An analytic approximation for the discrete frequencies 
cr„ can be found in [48]. We used this approximation as 
an initial input value for the frequency cro of the funda¬ 
mental mode n = 0, numerically integrating the radial 
perturbation equation (2.11) and searching for the value 
of CTo for which the scalar field perturbation vanishes on 
the mirror. 

Figure 1 shows the real and imaginary parts of the 
fundamental mode frequency ctq as a function of mirror 
radius for a selection of black hole and field charges, 
Q and q. The plots illustrate the following point: with 
the mirror placed close the black hole, the n = 0 mode 
decays exponentially; with the mirror placed far from the 
black hole, the n = 0 mode grows exponentially, generat¬ 
ing a superradiant instability; between these regimes is a 
‘transition point’, at exactly a = ac, at which the scalar 
field is in equilibrium with the black hole. 

In this paper we restrict our attention to a mass¬ 
less charged scalar field, but the superradiant insta¬ 
bility shown in Fig. 1 is also present when a massive 
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charged scalar field is considered [45-47]. In Ref. [45], a 
frequency-domain analysis was performed and a super¬ 
radiant instability found for massive charged scalar field 
modes with i = m = 1, where i is the total angular mo¬ 
mentum mode number and m the azimuthal mode num¬ 
ber. A time-domain study was undertaken in [46], again 
for a massive charged scalar field. In the £ = 1 case, 
the results of [46] show that at late times, the fundamen¬ 
tal unstable mode found in [45] dominates the evolution. 
They also find a superradiant instability for the £ = 0 
(spherically symmetric) mode, which grows more quickly 
than the £ = m = 1 unstable mode. The growth time of 
the £ = 0 unstable modes of the massless charged scalar 
field (whose frequencies are shown in Fig. 1) is of a sim¬ 
ilar order of magnitude to those in [46] for a charged 
scalar field with mass O.IM. The highly-explosive (yet 
still linear) regime was studied in Ref. [47]. 


C. Static hairy black holes 


At first glance, perturbations at the critical frequency 
tr = (Tc are time-dependent via ansatz (2.10), and thus 
not static. However, we may choose a gauge in which the 
scalar field is static, by inserting y = erf into Eq. (2.5). 
This gauge transformation removes the time-dependence 
from the field, and introduces a static constant term to 
Aq. This raises the possibility that static solutions may 
also exist for nontrivial scalar field $. 


1. Field equations 


To investigate this possibility, we now consider the 
spherically symmetric spacetime defined as follows 

ds"^ =-fhdt'^ + f-^dr"^ Fr'^dn'^, (2.13) 

where / = /(r) and h = h{r) and dFl^ is given by (2.8). 
We may write 


Hr) , 1 _ Mr) (2^14) 

r 

where m = m(r) is interpreted as the total mass within 
the given radius r. We assume that the static scalar field 
is real and depends only on the radial coordinate r, set¬ 
ting $ = (j}{r). Since we are considering a spherically 
symmetric spacetime, we can set the Ag and A^p com¬ 
ponents of the electromagnetic gauge potential to zero, 
and then use a gauge transformation (2.5) to set A^ to 
vanish. Thus the electromagnetic gauge potential takes 
the form A^ = [Ao(r), 0,0, 0]. 

With the above ansatz, the equations of motion (2.2) 


yield four non-trivial equations, 


h' = 
kE^ = 


0 = 

0 = 


(^) 


/\2 




f'h + -fh' -I-(/ — 1) 

2 r 


2/ fh' 


/A" + ^ ^ A' - q^^A 


2h 


2 / 


fh' 




{gApf 

fh 


(2.15a) 

(2.15b) 

(2.15c) 

(2.15d) 


where k = SttG and = (Ag)^. A prime ' denotes d/dr. 
For 0^0, Eqs. (2.15) cannot (apparently) be solved 
analytically, and some numerical analysis is required. 


2. Boundary conditions 


Let us now consider appropriate conditions to impose 
on the fields at the black hole horizon {r = r^) and at 
the mirror (r = r^). 

We assume that there is a regular event horizon defined 
by f{rh) = 0 and f'{rh) > 0. Thus, rrih = mfrh) = 
\rh. We demand that all physical quantities are finite 
in a future-horizon-penetrating coordinate system. This 
implies that the vector potential is zero at the horizon, 
A{rh) = 0; and h{rh) is finite. Without loss of generality, 
we set h{rh) = 1, which corresponds to a gauge choice in 
the definition of the time coordinate t. The scalar field 
equation (2.15d) implies that (f'frh) = 0. Hence regular 
Taylor series expansions of the field variables about r = 
r/i take the following form 

m = mh + m'^ir - rh) -f 0(r - ruf', 
h = \ + h'nir - rh) + 0{r - rh)"^, 

Ao = Eh{r - rh) + ^{r - rh)"^ + 0{r - rh)^, 

(l) = (l)h + ^{r-rh)'^+ 0{r-rh)^, (2.16) 


where Eh = Ag(r?t) is the electric field on the horizon. 
Inserting these expansions back into the field equations 
(2.15) gives 


m'h 

K 


A 


// 

h 




4 ’ 

Anq^rlcgEl 


^ r 2q\lcl>l 

[{KrlEl-2) 

{KrlEl-2)^- 


2 



(2.17) 


At fixed q,rh and Eh, these expansions are determined 
by just one further constant (ph- 
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FIG. 2. Left: an example plot of the two metric functions /, h and two matter functions Aq, rj) for a particular static black 
hole solution with q — 0.9, <f)h = 0.4 and Eh = 0.8. Right: Scalar field profiles for three different black hole solutions for fixed 
cfih = 0.3 and Eh ~ 0.6 and three values of the scalar-field charge q. 
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FIG. 3. An example plot of different black hole solutions with scalar charge q = 0.1. Left: Three scalar field profiles which 
share the same location of their first node at Vm ~ 27. Right: Three scalar field profiles with fixed Eh = 1 with a common 
node; their first (red, solid), second (blue, dashed) and third (green, dot-dashed) nodes coincide at ~ 27. 


We insist that the scalar field vanishes at the location 
of the mirror 

(l>{rm) = 0. (2.18) 

No further conditions are applied at r = r^, so the metric 
functions /, h and the electric gauge potential Aq are 
unconstrained at the mirror’s location. 


3. Solutions 

We seek static solutions by numerically integrating the 
equation set (2.15). To avoid the regular singular point at 
r = rt, we use the series expansions (2.16-2.17) as initial 
data, evaluated at r = r/j-be (where typically e ~ 10“^^). 
We choose units such that k = 8ttG = 1. 

Without loss of generality, one may rescale all dimen¬ 
sionful quantities by to obtain a dimensionless equa¬ 
tion set. Equivalently, one may simply set = 1 (thus 
mh = 1/2), leaving three free parameters: q, the scalar- 


field charge; (ph, the scalar-field magnitude on the hori¬ 
zon; and Eh = A'Q{rh), the electric field on the horizon. 
Henceforth, these should be thought of as dimensionless 
quantities, in units of Vh- We now explore this three- 
parameter solution space. 

The left-hand plot in Fig. 2 shows the four field vari¬ 
ables /(r), h{r), Ho(r) and (j)(r) for the case of scalar 
field charge q = 0.9, iph = 0.4 and Eh = 0.8. The scalar 
field (j){r) oscillates about zero; the other three field vari¬ 
ables f(r), h{r) and 2 lo(r) are monotonically increasing 
with r. Note that we do not expect that Ao{r), f{r) and 
h{r) will necessarily approach finite limits as r —>■ oo, 
since there are no asymptotically flat black hole solu¬ 
tions in this model with nontrivial scalar field hair [31]. 
The right-hand plot in Fig. 2 shows an example of the 
scalar field profile outside the horizon for three values of 
q with fixed 4>h =0.3 and Eh =0.6. Here, the oscillating 
behaviour of the scalar field is more clearly seen. One 
could obtain a black-hole-in-a-cavity solution by placing 
the mirror at any of the nodes of the scalar field. For the 
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FIG. 4. Numerical exploration of the parameter space for static solutions with scalar field hair in a cavity. The system has 
three free parameters: q (the charge of the field), (ph (the scalar field on the horizon), and Eh (the electric field on the horizon). 
The plots show two-dimensional slices of the solution space with (top row, left to right) g = 0.1, q = 0.2, (bottom row, left to 
right) q = 0.4 and q = 0.8. The shaded area indicates the region where solutions exist, with the scalar field having at least one 
node, and /(r) > 0 for r > rn- Note that solutions exist throughout the central region all the way towards (j)h —t 0, Eh —>■ 0; 
except along the line Eh = 0. The coloured lines are contours of constant mirror radius rm, where Vm lies at the first node of 
the scalar field (j>. 


majority of this paper we will consider the case that the 
mirror is located at the first node. 

It is possible to have different black hole configura¬ 
tions which share the same mirror radius, as illustrated 
by Fig. 3. In the left-hand plot in Fig. 3, we show three 
scalar field profiles (f>{r) for different values of the elec¬ 
tric field at the horizon Eh and fixed scalar-field charge 
q = 0.1. Each solution, despite having different values of 
(ph and Eh, has a node of the scalar field at ~ 27. A 
further three distinct scalar field profiles are displayed in 
the right-hand plot in Fig. 3. For a given Eh = 1 (the 
scalar-field charge is still fixed to be g = 0.1), the first, 
second and third nodes of the scalar field share the same 
location. 

Figure 4 illustrates the three-dimensional solution 
space of static hairy black holes in a cavity. The plots 


indicate that solutions exist in a contiguous region of 
{q, Eh, (ph} parameter space, in which solutions with at 
least one node in the scalar field are permitted. Out¬ 
side this region, we find that an excess of stress energy 
causes the metric function /(r) to develop an additional 
zero before the scalar field develops its first node (sug¬ 
gesting that an additional horizon forms). We note that 
(i) solutions with <ph = 0 are well-known: these are the 
Reissner-Nordstrom black holes; and (ii) uncharged so¬ 
lutions, with Eh = 0 or q = 0, are not possible, as the 
scalar field does not develop a node. 
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III. STABILITY ANALYSIS 

In the previous section, we learned that a fully coupled 
system of gravity and a charged scalar field admits black 
hole solutions with scalar hair confined within a cavity. 
Are such hairy configurations stable or unstable? If such 
solutions are stable, it is plausible that they represent the 
end-point of the superradiant instability for a massless 
charged scalar perturbation on the Reissner-Nordstrom 
background with a mirror, as described in Sec. IIB. In 
this section we examine the stability of the hairy black 
hole configurations under linear, spherically symmetric, 
perturbations. 


A. Dynamical equations 


We begin by outlining the ansatz for the field variables. 
We consider a spherically symmetric metric of the form 
(2.13), but with the functions f{t,r) and h{t,r) depend¬ 
ing on time t as well as the radial coordinate r. The 
scalar field d) = (j>{t, r) similarly is time-dependent but is 
now complex. By virtue of spherical symmetry, we can 
set Ag = 0 = A^. Using (2.5), we may apply a gauge 
transformation to eliminate leaving A^ with only a 
temporal component, = [Ao(t, r), 0, 0,0]. 

It is convenient to introduce a new metric variable, 

7 = (3.1) 


From the Einstein field equations (2.2a), and, in par¬ 
ticular, from the combinations Gtt, Gu + J^Grr and 
Gtt — j'^Grr-, we obtain 


11 . 

f 

h' 

ll 

7 


rr 

7 

+ (!-/)> 


(3.2a) 

(3.2b) 

(3.2c) 


where 


r = -b |7^'P + + 2gAoIm((/)0*). (3.3) 

Here, the dot ’ and prime ' denote partial derivatives with 
respect to t and r, respectively, and the asterisk * denotes 
complex conjugation. We note the equation for 7 has no 
dependence on r, and thus it does not explicitly depend 
on the scalar field </>. 

There is one further independent, nontrivial compo¬ 
nent of the Einstein field equations (2.2a), namely the 
Gtr component, which gives 

-j= rRe + rgAoIm . (3.4) 


Note that in the static limit this component is identically 
zero. 

From the t and r components of the Maxwell equations 
( 2 . 2 b) we obtain two dynamical equations 

= Jt = Ao - qlui , (3.5a) 

1 /rA' \ 

{jA^J " = -9lm(70'</>*)- (3-5b) 


The scalar field equation (2.2c), first written in the 
form rfhO’^Dacj) = 0 , yields 


0 = — 4) + -Ip + 7 — ■^tp + 2iqAoip + iqAopj 

7 r 


- iq-Aoip + q^Alp), 


(3.6) 


where 


Ip = rep. 


(3.7) 


The equations (3.2-3. 6 ) govern how the spacetime met¬ 
ric, electromagnetic held and the massless scalar held 
evolve with time. 


B. Perturbation equations 


Our aim in this section is to study the stability of the 
hairy black holes found previously in Sec. IIC. We there¬ 
fore now consider linear perturbations around a non¬ 
vacuum solution by introducing the following notation: 


f = f{r) + Sf{t,r), 
h = h{r) -b 5h{t, r), 

7 = i{r) +6j{t,r), 

^0 = ^o(’') + SAQ{t,r), 

Ip = Tpir) + Stp{t, r). (3.8) 


In this formalism, / is the equilibrium quantity and Sf 
is the perturbation. We assume that only 6ip is a com¬ 
plex variable while all other quantities are real. From 
(3.2, 3.4-3. 6 ) six independent dynamical equations can 
be obtained. For the remainder of this section, we work 
to hrst order in the perturbations. 




(3.9a) 


First, (3.2a) gives 
Sf 


1 , rjA'^r r 

rp 272 / 2 _ 


sf = - 




' K A' _ 

v2 


iqAotp 


-2 *^210 + 

'P'y'^ /yU 


q^AlP 


rf{A',f 


^7 


I - sp'J 


+ 


2rp 

ii)' i’q^Al {[} 

2r^ 


Ip tp' 


{Sp + Sp*) 


2 r 2 2rp 

A similar equation for Sh can be obtained from (3.2b), that is, 


(S-ip + Sp). 


Sh' h' 2AQq^p ^ ^ ^ 2q^AQp ^ pA^ip 


x-f“- 


rp 


-SAo 


rpryO 


-Sj = 




{si; - Sp'^ 

■0 Tpq^AQ p)' 


rp 


t-i} (iv' + v) 

(^0 + Sip*). 


(3.9b) 


It will be useful for our later analysis to have an equation with the same structure for S^. This can be derived from 
(3.2c), leading to 


Sp 


1 


7 P 


rim 


■7 


rJK 


57 + 


1 


rA' 




(3.9c) 


Note that Eq. (3.9c) is not an independent equation because it can be derived directly from the definition of S^ = 
S{fP^^) and Eqs. (3.9a-3.9b). The final component of the Einstein field equations (3.4) takes the form 


SJ 

f 


Ip' Ip 

2r 2P 


l + sp'^ 




2r 2r 

The two components of the Maxwell equations (3.5) yield the following expressions, 


(3.9d) 


y/h 


<5A" + 


y/h 


2 h' 

r 2h 


. . y , , 7A(, , 1 

SAq — —07 "I—- f=Sh + ^ 

P Pi 2hVh 2h 


iP hy/h 


Sh 


SA', 


y/h 2hy/h 

Lastly, the Klein-Gordon equation (3.6) yields 


+ '£ip- p + (‘'ll + 

'’i sh=lll - sr) HI . 

2 r 2 


2P 


(3.9e) 

(3.9f) 


0 = —Sip + P Sip" + 2iqAoSip + 'ypSip' + 

07' 


„2 12 77 

1 ^ 0 - — 


Sip — 


iqApip 


Sj ■ 


70' - 


70 


Sj' 


200 " + pip' - 


^7 + iqipSAp + 2q ipApSAp. 


(3.9g) 


It can be seen from (3.9d, 3.9f) that the imaginary part of the scalar field perturbation, Im((50), is out of phase 
with Sf,Sh,SAo and the real part of Sip. Eor this reason, we decompose the perturbed scalar field in the following 
way 


Sip(t, r) = Su(t, r) + iSw{t, r), 


(3.10) 


where u{t,r) and w(t,r) are real perturbations. This definition implies that Sw(t,r) is only determined up to an 
arbitrary function of r. With the definition (3.10), the Klein-Gordon equation (3.9g) can be now separated into two 
independent equations, corresponding to the real and imaginary parts. 

By integrating once with respect to time t, (3.9d, 3.9f) imply that 


^ = 1 
/ r 

Sh 


0 7,' 

- iP 

r 




(3.IIa) 
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(3.IIb) 
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where 5!F{r) and 5'H{r) are arbitrary functions of the radial coordinate r. Thus it is straightforward to obtain 


Jy = 


7 






[f^ + + -^17 Sw' + -^6 Aq + + ^-4^8%. 


qii’ [IT 


7 


fh, 


r'^AL 


r'^A'r, 


A' 


(3.11c) 


Equations (3.11) allow us to rewrite the metric perturbations Sf,Sh and dy in terms of the matter perturbations 
Su,Sw and SAq, and hence eliminate the metric perturbations from the Maxwell (3.9e-3.9f) and Klein-Gordon (3.9g) 
equations. Moreover, it is possible to construct the following linear first-order differential equation, from (3.9a, 3.9e) 


5T' 


The above equation is integrable, with solution 


If h' 1 

— -b ^ -b - 


'q^oVhf;^ ^A'4 fA'oh'' 

[/ 2h rj 

2y 2y2 

Ao 2'/h 


sn. 


5F = 

27 


(3.12) 


(3.13) 


up to an overall constant. We can use this relation to eliminate 6% from the perturbed field equations. 

To find the equations for the matter perturbations, we begin with (3.9a), obtaining the following equation once the 
metric perturbations have been eliminated: 


0 = dui — "f‘^6w" + 

+ qAo 


- 77 ' ' 




2 A' 


A 


_9 _L ^ _ tltl. 
^ ' 2 

r 


Su - 


Sw' + 

qAo'4 


„2 y2 , fh 7^ fA'4 q^AoiPtlj' ^ 

q -^0 ' 0 9 n n * / 


271' 


6w 


A' 


SAq- qipSAo 


ATT' , 

-T— j8~r -b 

qAoip 


fh 


fAA'4 , 


qAotjj 2qAo4i qipAl rA[j 




(3.14) 


where 


A = fh -b rA qA'q. 

The imaginary part of scalar field equation (3.9g) can be integrated once with respect to time to give 

q^Aofff’' 


(3.15) 


0 = (5?h — 'y Sw + 


__ , 
-77 H- ^ T, A 


AA', 


+ qAo 


-2 


Tp'^ 'iplp' 


Su - 


Sw' + 
qAofi 


-n^ - 




-A 


li' 


Sw 




SA'q - qipSAo + sg{r), 


(3.16) 


where SQ{r) is an arbitrary function of the radial coordinate r. We next compare the two equations (3.14, 3.16). This 
gives another linear first-order equation, this time for S!F and SQ: 


Q = ST' + 


A', A[ 


10 


/ 


ST- 


qAoip 

ry^ 


sg. 


(3.17) 


We will return to Eq. (3.17) in the next section where we eliminate the unknown functions ST{r) and Sg{r). For the 
last step in our derivation of the linearised perturbation equations, we use (3.9b, 3.9c) to eliminate Sw and SAq from 
the real part of the Klein-Gordon equation (3.9g). 

Following these steps, we may obtain three equations governing three perturbations: Su, Sw and SAq. The first 
equation is derived from the real part of the Klein-Gordon equation (3.9g) and takes the form 


0 = Sil — j'^Su" — yy'du' -b 


3q‘^Af 


77 ' 


-fh 


Tx '2 


+ qf Aq 

2 ry 2 

-I- ^ST -b y 

if 


2^/^y+G4!^A + ^ 


' 


h rA'^ 




fA'4 


Sw' + qAo 




277' 


2q^Al-^+j4' 


Tin 


-/ 7 

■7 - 

r 


Sw 


I 


in)' - 


nA'o 


ST. 


(3.18a) 
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The second equation is obtained from (3.16), after application of (3.17) 

q'^ Ao’tjj'ij}' 


0 = 6w — "/^Sw" + 
, 9^0 '0 


__ ■ 


5w' 


-0^4^ - 


r^A'n 


-A+ — 


Sw — qAo 


2 + r\* 


A' 


SA',-qi^SAo-^S:F'' 


qAat/j 


qAotjj 


1 r A' A'' 

---T + ^ + ^-r 


f bio 


bl' 


Ts / 2 ' 

il; 


ST. 


Su 


(3.18b) 


The third equation comes from the Einstein field equation (3.9c) 


0 = AKaSw" + Ttr 


A'or^ 




V ' 


Su' - 


2lo^o7 




r^A';^ 


A 


2 rAoA',j ( ^ + 


rg 2 ^ 2 \ q^h'il;'ip' 


r'^A'i^ 


Sw 


1 +b 

r 7 




Su + 


^bl' 


bl; 


0 J 


A 


rAoA'n 


ST 


rf Ao 2fhj 


To summarise, we have obtained two dynamical equa¬ 
tions (3.18a, 3.18b) which involve time derivatives, and 
one constraint equation (3.18c) which contains only 
derivatives with respect to r. Note that that all these 
equations (3.18) only contain radial derivatives of the 
electric potential perturbation SAq, and not time deriva¬ 
tives. Essentially, this is due to residual gauge freedom 
(as discussed in Sec. Ill A), which means that an arbi¬ 
trary global function of time can be added to SAq without 
changing physical quantities such as the electromagnetic 
field. 


C. Boundary conditions 


We now consider the boundary conditions for the per¬ 
turbed field variables at two boundaries, i.e. at the black 
hole horizon (r = rh) and at the mirror (r = r^). Near 
the horizon, we impose ingoing boundary conditions 


Su{t, r) = Re 
Sw{t, r) = Re 
^Ao(t, r) = Re 


e-*'"(‘+’’*)i(}(r) 

e-*'"(‘+’’*)io(r) 


where r* is the tortoise coordinate defined by 
dr^. 


dr 


= 7- 


(3.19) 


(3.20) 


Here, ft, w and Aq are complex functions which depend 
only on the radial coordinate r and have Taylor series 
expansions near the horizon of the form 

M = Mo -I- Mi(r - rh) -I- 0{r - r/^)^, 
w = Wo + wi{r - rh) + 0{r - rh)^, 

Aq = Ai{r- rh) + A 2 {r - r/,)^ -f 0{r - rh)^. (3.21) 


Before we proceed further, we noted earlier that adding 
an arbitrary function of r to Sw makes no difference to 


the scalar field perturbation Sip- This freedom allows us 
to set (55 = 0 in (3.16). Hence, (3.17) is solvable using the 
conventional integrating factor method and the solution 
is given by 


ST = IC 


rApA'^ 

f 


exp 



(3.22) 


where /C is a constant of integration. As ingoing bound¬ 
ary conditions are required for all perturbations, includ¬ 
ing the metric variables Sf, Sh (and thus ( 57 ), it must be 
the case from Eq. (3.11a) that ST = 0 at r = rh- There¬ 
fore we must set /C = 0 in (3.22), and so ST vanishes 
identically. Thus ST' and ST are eliminated from the 
perturbation equations (3.18) by our choice of boundary 
conditions. 

By inserting (3.19, 3.21) into the perturbation equa¬ 
tions (3.18), we find that ui, wi, Ai and A 2 can be ex¬ 
pressed in terms of up, wp and a. The simpler of these 
expressions are 

7 iga(/>h (2 -f rlEf) 

^ rh{-2 + 2irh<7+ rlEl)^^' 

~ _ -SqEha^rlwp + (-2 -f rlElf up 
^ rh (-2 -f rlEl) (-2 -f AirhO + rlE^) ’ 

while the expressions for wi and A 2 are sufficiently com¬ 
plicated to be omitted here. Thus with given values for 
the background parameters q^fph a-nd Eh, the boundary 
conditions (3.19) depend on three additional parameters, 
namely up, wp and cr. We emphasize that the parameters 
itp, Wp and a are all complex. The physical perturbations 
arise from taking the real part in (3.19). 

At the mirror r = r^, the scalar field perturbation Sip 
(like the background scalar field ip) must vanish. The 
perturbations of the metric functions and electric poten¬ 
tial are unconstrained there. Since the real and imagi¬ 
nary parts of the scalar field perturbation (3.10) take the 
form (3.19), at the mirror the functions u(r) and w{r) 
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FIG. 5. An example plot of the three perturbation functions, u, w and Aq satisfying the perturbation equations (3.18) 
and boundary conditions (3.19, 3.24). The equilibrium solution parameters are: scalar-field charge q = 0.1, electric field on 
the horizon Eh = 0.8 and scalar field on the horizon (ph = 1.2. The corresponding eigenvalue a is 0.1731 — 0.0038i and 
uq = 0.4397 -I- 0.02317 The horizontal and vertical lines are included to help visualise where the perturbation functions vanish 
and the location of the mirror. For this example the mirror is situated at Tm ~ 24, which is the first zero of the equilibrium 
scalar field. 


must satisfy 

u{rm) = Q = w{rm). (3.24) 

We require both the real and imaginary parts of u{r) 
and w{r) to vanish at the mirror so that the real and 
imaginary parts of the scalar field perturbation (given 
by (3.19)) vanish for all time t, when the real part in 
(3.19) is taken. 

In summary, using the form (3.19) for the matter per¬ 
turbations, we now have three ordinary differential equa¬ 
tions (3.18) for three unknown functions of r, namely: 
u, w and Aq. Together with the boundary conditions 
(3.19, 3.24), we now have a system which can be solved 
numerically. 

D. Method and results 

We implement a shooting method to numerically solve 
the boundary value problem (3.18, 3.19, 3.24). Since 
both the perturbation equations and boundary condi¬ 
tions are linear, we set the overall scale of the pertur¬ 
bations so that wq is fixed to be unity. This leaves two 
free parameters, uq and a, which we use as shooting pa¬ 
rameters. The process of numerical integration is as fol¬ 
lows. Firstly, we specify the background parameters q, (ph 
and Eh, then integrate the static field equations. We ob¬ 
tain the numerical hairy black hole solution and find the 
location of the first zero of the equilibrium scalar field, 
setting this to be the mirror location Vm- Secondly, the 
three coupled perturbation equations (3.18) are solved by 
searching for values of uq and a such that the boundary 
conditions (3.19, 3.24) are satisfied. We are particularly 
interested in the sign of the imaginary part of the fre¬ 
quency, Im(tT). Perturbations for which Im(CT) < 0 are 
stable and decay exponentially in time, whereas pertur¬ 
bations for which Im((T) > 0 are unstable, growing expo¬ 
nentially in time. 

As an example, we plot in Fig. 5 the behaviour of 
u, w and Aq for scalar charge q = 0.1, with horizon 


values for the electric field Eh = 0.8 and scalar field 
(ph = 1.2. The values of the two shooting parameters 
are a = 0.1731 — 0.00381 and uq = 0.4397-1-0.02311. This 
figure clearly demonstrates that both the real and imag¬ 
inary parts of the field variables u and w vanish at the 
location of the mirror. By contrast, the perturbation of 
the electric potential Aq does not vanish on the mirror. 
This particular perturbation mode decays exponentially 
with time since the frequency a satisfies Im(cr) < 0. 

The key question we explore in this section is whether 
the perturbation mode shown in Fig. 5 is typical in having 
Im(o-) < 0. As illustrated by Fig. 4, we have a three- 
dimensional parameter space of equilibrium solutions, 
governed by the parameters q (the scalar-field charge), 
(ph (the value of the scalar field on the horizon) and Eh 
(the electric field on the horizon). It is clearly impracti¬ 
cal to test every possible solution in this phase space, and 
therefore in this section we present a selection of results 
probing various parts of the phase space. 

With the mirror at the first zero of the equilibrium 
scalar field, for fixed values of the parameters q, (ph and 
Eh, we search for perturbations solving the equations 
(3.18), subject to the boundary conditions (3.19, 3.24). 
For each fixed equilibrium solution we found a single 
value of a (together with one value of the other shooting 
parameter uq) such that the corresponding solution of 
the perturbation equations satisfies the boundary condi¬ 
tions. In all cases examined with the mirror at the first 
node, we found that Im((T) < 0, implying that the per¬ 
turbations exponentially decay in time, suggesting the 
equilibrium solutions are stable. 

We now present a selection of numerical results. In 
Figs. 6-10 we fix two of the parameters q, (ph and Eh 
and vary the third. 

Figure 6 shows the real (left) and imaginary (right) 
parts of the frequency cr as a function of mirror radius 
Tm (top row) and the value of the equilibrium scalar field 
(ph (bottom row) for fixed g = 0.1. In each plot, three dif¬ 
ferent curves represent three distinct values of Eh, and 
on each curve (ph varies from 0.1 to 1.4. Our numeri- 
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FIG. 6. The real (left) and imaginary (right) part of the mode frequency cr is plotted as a function of (top row) the mirror 
radius r-m and (bottom row) the equilibrium scalar field on the horizon (j>h, for q = 0.1, 4>h G (0.1, 1.4) and various values of 
Eh- The mirror is located at the first zero of the equilibrium scalar field. In all these plots, Im((T) < 0, and the perturbations 
decay exponentially in time. 


Di 




El, Eh 


FIG. 7. The real (left) and imaginary (right) part of the mode frequency a is plotted as a function of the electric field at the 
horizon Eh, for g = 0.1 and various values of (j>h- The mirror is located at the first zero of the equilibrium scalar field. In the 
right-hand plot, lm{a) < 0 and we have plotted the logarithm of the modulus of Im((T). 


cal method breaks down when cph is very small and the 
mirror location is far from the black hole event horizon. 
When (j)h is greater than ~ 1.4, no static black hole solu¬ 
tions exist with nontrivial scalar field hair (see Fig. 4); in 
the case where Eh = 0.8, static hairy black holes cannot 
be found for (ph larger than ~ 1.2. As can be clearly seen 
in the Im(CT) plots, all black hole solutions examined ap¬ 


pear to be linearly stable against spherically symmetric 
perturbations. 

Figure 6 shows how Im(CT) decreases (so that the per¬ 
turbations decay more rapidly) and Re(CT) increases as 
(fih is increased, and the mirror moves closer to the black 
hole horizon. For the values of Eh shown in this figure, 
Im((T) also decreases as Eh increases for fixed (ph- 
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FIG. 8. The real (left) and imaginary (right) part of the mode frequency a is plotted as a function of the scalar charge q, for 
various values of (ph and Eh- The mirror is located at the first zero of the equilibrium scalar field. All values of Im((T) shown 
in the right-hand plot are negative. 
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FIG. 9. The real (left) and imaginary (right) part of the mode frequency a is plotted as a function of the electric field at the 
horizon Eh, for fixed (j)h = 0.2 and various values of q. The mirror is located at the first zero of the equilibrium scalar field. 
All values of Im((T) shown in the right-hand plot are negative; we have plotted the logarithm of the modulus of Im(cr). 


Figure 7 illustrates the real (left) and imaginary (right) 
part of the frequency cr against the electric field at the 
horizon Eh for three distinct values of (j)h- It remains 
the case that Im((T) < 0; as it can take a small value, 
we plot the logarithm of the modulus of Im(cr). As the 
electric field on the horizon, Eh, increases for fixed tph, 
then Re(CT) increases and Im(CT) decreases. Furthermore, 
as (j)h increases, Re(CT) increases and Im(tT) decreases. 

Let us now explore the effect of changing the scalar- 
field charge q. We may fix the values of iph and Eh, and 
vary q. Figure 8 shows Re((T) (left) and Im((T) (right) as 
functions of the charge of the scalar field q. Once again, 
we find only stable modes with Im(cr) < 0. As the scalar- 
field charge q increases, we see that Re(cr) increases and 
Im((T) decreases. 

In Fig. 9, four different scalar-field charges are chosen, 
and the value of the equilibrium scalar field on the hori¬ 
zon is fixed to be (j)h = 0.2, with Eh varying. Notice that 
as q increases the real part of the frequency a increases 
and the imaginary part decreases (we plot the magnitude 
of Im((7) on a logarithmic scale as it is small). Moreover, 


from Fig. 9, we observe that as q decreases, the curves 
for Re((T) and Im((7) cover a greater range of Eh values. 
This is due to the fact that as q becomes smaller the two- 
dimensional phase space of static solutions expands (see 
Fig. 4). 

Finally, in Fig. 10 we display the real (left) and imag¬ 
inary (right) part of the frequency cr as a function of 
mirror radius (top row) and the value of the equi¬ 
librium scalar field on the horizon (f>h (bottom row) for 
fixed scalar charge q = 0.2. This should be compared 
with the corresponding plot in Fig. 6 for g = 0.1. We see 
in Fig. 10 that as Eh increases the range of values of (ph 
also increases. This can be understood from the phase 
space plot of the static solutions (upper-right-hand plot 
in Fig. 4). Again, Im(cr) remains negative. As in previ¬ 
ous figures, with fixed Eh, increasing (ph increases Re(cr) 
and decreases Im(cr), while increasing Eh for fixed pt 
decreases Im(cr). 

In all the figures considered so far in this section, the 
mirror was located at the first zero of the equilibrium 
scalar field. We have found a consistent picture: for each 









14 






FIG. 10. The real (left) and imaginary (right) part of the mode frequency a is plotted as a function of (top row) the mirror 
radius Vm and (bottom row) the equilibrium scalar field on the horizon 4>h, for q — 0.2, <f)h £ (0.1,1.3) and various values of 
Eh- The mirror is located at the first zero of the equilibrium scalar field. In the right-hand plots, we plot the logarithm of the 
modulus of Im((j); note however that all values of Im(CT) shown are negative. 


static hairy black hole, we can only find perturbations 
which decay exponentially in time. We conclude that 
the static hairy black holes with the mirror at the first 
zero of the equilibrium scalar field appear to be stable. 

We close this section by considering some results when 
the mirror is located at the second zero of the equilibrium 
scalar field, providing an example plot in Fig. 11. The 
real (left) and imaginary (right) part of the frequency cr 
are plotted against mirror radius (top row) and the 
value of the scalar field on the horizon (j)h (bottom row) 
for fixed g = 0.1. In contrast to the hrst-zero case, we 
now find perturbations with Im(cr) > 0 for all the static 
black hole solutions considered in Fig. 11, so that the 
perturbations grow exponentially in time. We conclude 
that static hairy black holes with the mirror at the sec¬ 
ond zero of the equilibrium scalar field are unstable. We 
conjecture that if the mirror was located at a node after 
the second zero of the equilibrium scalar field, then the 
hairy black holes would remain unstable. 


IV. CONCLUSIONS 

When a charged scalar field interacts with an electri¬ 
cally charged Reissner-Nordstrom black hole surrounded 


by a reflecting mirror, superradiantly-unstable modes ex¬ 
ist [45, 46]. A natural question arises: what is the end¬ 
point of this instability? This question has been the focus 
of our work in this paper. 

Working in the frequency domain, we have confirmed 
the time-domain results of [46], namely that the charged 
scalar field, linearized around the $ = 0 background, has 
spherically symmetric unstable modes, if the mirror is 
sufficiently far from the black hole horizon. This led us 
to consider non-linear spherically symmetric solutions of 
the fully coupled Einstein-charged scalar field theory as 
possible end-points of this superradiant instability. In the 
‘scalar electrodynamics’ model, a charged complex scalar 
field is coupled to an electromagnetic field and the usual 
Einstein-Hilbert gravitational Lagrangian. Solving the 
equilibrium field equations, we found black hole solutions 
with a nontrivial scalar field which oscillates about zero. 
We may place a reflecting mirror at any one of the nodes 
of the scalar field, to obtain a black hole in a cavity. It 
is important to note that the solutions we find do not 
contradict the no-hair theorem of Bekenstein [31] which 
applies in the absence of a mirror-like boundary. 

To investigate whether these hairy black holes could 
be possible end-points of the superradiant instability, we 
have considered spherically symmetric perturbations of 
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FIG. 11. Instability of perturbations when the mirror is placed at the second zero of the static scalar field. The real (left) 
and imaginary (right) part of the mode frequency o is plotted as a function of (top row) the mirror radius Cm and (bottom 
row) the equilibrium scalar field on the horizon ij>h, for q = 0.1, (j)h G (0.1,1.4) and various values of Eh- Note that in these 
plots, unlike the plots for the first-zero case, the imaginary part of the frequency is positive, implying that the perturbations 
are exponentially growing in time. 


the hairy black hole solutions. These perturbations sat¬ 
isfy ingoing boundary conditions on the horizon. Fur¬ 
thermore, the perturbations of the charged scalar field 
vanish on the mirror. The resulting perturbation equa¬ 
tions, though linear, are highly coupled and can (we be¬ 
lieve) only be solved numerically. 

With the mirror placed at the first zero of the equilib¬ 
rium scalar field, we find no evidence of any instabilities: 
all such equilibrium solutions appear to be stable under 
small perturbations. On the other hand, if the mirror is 
placed at the second zero of the equilibrium scalar field, 
we find unstable perturbations which grow with time. 
We conjecture that the same result would hold if the 
mirror were situated at the third or subsequent zeros of 
the equilibrium scalar field. 

Let us now consider the implications of these results 
in a wider context. It is known that a superradiant in¬ 
stability can also arise when a charged black hole is em¬ 
bedded in a asymptotically Anti-de Sitter (AdS) space- 
time [48, 49]. In such a scenario, the timelike boundary 
of spacetime itself acts as the mirror, providing a natu¬ 
ral reflecting boundary condition (or ‘Dirichlet wall’). In 
Refs. [50, 51] asymptotically AdS ‘hairy’ black holes were 
constructed within the context of supergravity/higher¬ 


dimensional theories. Two influential ideas underpin¬ 
ning this work can be traced through Refs. [41, 48-51]: 
(i) ‘hairy’ stationary solutions are plausible endpoints 
of the superradiant instability; and (ii) in the limit of 
small field amplitude, the hairy solutions connect to (per¬ 
turbed) vacuum solutions endowed with linear perturba¬ 
tions at the critical superradiant frequency It is plau¬ 
sible that the stability properties of the four-dimensional 
hairy black holes, considered here, may be shared by 
their higher dimensional, non-asymptotically flat cousins. 
However, this remains to be investigated. 

Our view is that the ‘scalar electrodynamics’ model 
may also provide insight into the development of super¬ 
radiant instabilities in astrophysical systems. However, 
we should proceed with an element of caution, for two 
reasons. First, in the Kerr black hole case, superradi¬ 
ance is promoted by angular momentum, rather than 
charge. Thus, the Kerr instability does not appear in 
the spherically-symmetric sector, making any stability 
analysis considerably more involved. Second, in the Kerr 
case, instabilities can arise spontaneously in bound states 
of an (ultra-light) massive bosonic field. By contrast, in 
the charged case, the competition between gravitational 
attraction and electrostatic repulsion means that bound 











16 


states cannot form in the superradiant frequency regime; 
thus the artifice of a mirror is necessary. These factors 
suggest that one should be cautious when attempting to 
infer from analogy. 

This work was motivated, in part, by the recent dis¬ 
covery of a Kerr-scalar family of (asymptotically-flat and 
four-dimensional) black hole solutions [39, 40]. The dis¬ 
covery inspired Herdeiro and Radu [52] to propose the 
following conjecture: ‘a (hairless) black hole which is 
afflicted by the superradiant instability of a given field 
must allow hairy generalizations with that held’. (As 
noted in Ref. [40], the given held should also generate a 
stress-energy with appropriate Killing symmetries; this 
excludes, for example, a real scalar held). Our equilib¬ 
rium charged-scalar black hole solutions (Sec. IIC) are 
in accord with Herdeiro and Radu’s conjecture. One 
could also envisage a stronger conjecture: ‘a (hairless) 
black hole solution afflicted by a superradiant instabil¬ 
ity of a given held may naturally evolve towards a hairy 
black hole solution which is stable under perturbations in 
that held’. When the mirror is located at the hrst zero 
of the equilibrium scalar held, our charged-scalar black 
hole solutions are stable (at least to time-periodic, linear, 
spherically symmetric perturbations). We have there¬ 
fore conjectured that they are possible end-points of the 
superradiant instability discovered in [45, 46], in accord 
with our stronger version of the Herdeiro/Radu conjec¬ 
ture. We should mention that the stronger conjecture is 
somewhat provocative, as it is apparently in tension with 
numerical simulations in the Kerr case which suggest that 
nonlinear effects lead to collapse of the held, and a subse¬ 
quent explosive phenomenon, known as a ‘bosenova’ [53] 
(for other perspectives see Refs. [54-58]). 

In order to ascertain whether these hairy charged- 
scalar black holes are indeed (i) stable under more generic 


perturbations, and (ii) natural end-points of the vacuum 
superradiant instability, a full nonlinear time-domain 
numerical simulation would be required. One could 
start with a small perturbation of a Reissner-Nordstrom 
black hole in a cavity, and track the development of the 
instability into the nonlinear regime. Such a nonlinear 
simulation would be a technical achievement, as (i) the 
horizon would be dynamical, and (ii) the growth rate 
of the superradiantly unstable modes of the Reissner- 
Nordstrom black hole in a cavity is about two orders of 
magnitude smaller than their frequency. Here, the meth¬ 
ods of numerical relativity may find a further application 
[59]. Such simulations could be greatly simplihed by 
restricting to spherical symmetry, which is not possible 
in the Kerr context. Alternatively, as a starting point, it 
would be instructive to perform a time-domain analysis 
of linear, spherically symmetric perturbations of the 
equilibrium charged-scalar hairy black holes presented 
in this paper, in order to confirm the frequency-domain 
stability results presented here. Either approach would 
surely lead towards a fuller understanding of the generic 
features of superradiant instabilities, and their relevance 
(or otherwise) in astrophysics. 
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